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Abstract 

The cross section of the diffractive process e + p —>• e + Xp is measured at a centre-of- 
mass energy of 318 GeV, where the system X contains at least two jets and the leading final 
state proton p is detected in the HI Very Forward Proton Spectrometer. The measurement is 
performed in photoproduction with photon virtualities Q 2 < 2 GeV 2 and in deep-inelastic 
scattering with 4GeV 2 < Q 2 < 80 GeV 2 . The results are compared to next-to-leading 
order QCD calculations based on diffractive parton distribution functions as extracted from 
measurements of inclusive cross sections in diffractive deep-inelastic scattering. 
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1 Introduction 


Diffractive processes, ep —> eXY, where the systems X and Y are separated in rapidity, have 
been studied extensively at the electron-proton collider HERA. In diffractive processes the in¬ 
teracting hadrons remain intact or dissociate into low mass hadronic systems via an exchange 
which has vacuum quantum numbers, often referred to as a pomeron (IP). Experimentally, 
diffractive events may be selected either by the presence of a large rapidity gap (LRG) in the ra¬ 
pidity distribution of the outgoing hadrons or by detecting a leading proton in the final state. The 
HI experiment was equipped with two dedicated detectors, the Forward Proton Spectrometer 
(FPS) 111] and the Very Forward Proton Spectrometer (VFPS) [0] to detect the leading protons. 

In the framework of the collinear factorisation theorem [0 diffractive parton distribution 
functions (DPDFs) may be defined. The factorisation theorem predicts that the cross section 
can be expressed as the convolution of non-perturbative DPDFs and partonic cross sections of 
the hard sub-process, calculable within perturbative Quantum Chromodynamics (QCD). The 
DPDFs have properties similar to the parton distribution functions of the proton, but with the 
constraint of a leading proton or its low mass excitations being present in the final state. 

DPDFs were obtained at HERA from inclusive diffractive deep-inelastic scattering (DDIS) 
data 0 , 0 . Given the DPDFs, perturbative QCD calculations are expected to be applicable 
to other processes such as jet and heavy quark production in DDIS at HERA Ml. Indeed, 
next-to-leading order (NFO) QCD predictions using DPDFs describe these measurements well. 


In diffractive hadron-hadron interactions however, the production of jets is found to be 
suppressed by about one order of magnitude HQ, as compared to predictions based on 
HERA DPDFs. This ’’factorisation breaking” may be explained e.g. by soft interactions or 
multi-pomeron exchanges between the hadrons and/or rescattering phenomena which destroy 
the diffractive event signature BOS. 


The issues of DPDF applicability and factorisation breaking can also be studied in hard 
diffractive photoproduction ( 7 p), where the virtuality of the exchanged photon Q 2 is close to 
zero. In the photoproduction regime, within the leading order approach, the small photon virtu¬ 
ality allows for partonic fluctuations 7 —> qq that last long enough to interact with the partons in 
the proton. In this regime the photon can be treated as a quasi-real target and therefore exhibits 
hadronic structure. 


Diffractive photoproduction of dijets in ep collisions at HERA have been measured by HI 
HE! and ZEUS [10. In each of these measurements diffractive events are selected by 
requiring a large rapidity gap. Different ratios of data to the NFO QCD prediction have been 
reported by HI and ZEUS: while HI reported their data to be suppressed by a factor of 0.6 with 
respect to the NFO QCD predictions II171 18], the ZEUS data are compatible with the theoretical 
expectations [20]]. Various mechanisms of suppressing diffractive dijet photoproduction have 
been proposed [ 21 , 22 ]. 

Enhanced sensitivity to the differences between theory and data may be achieved by calcu¬ 
lating the double ratio of the ratio of data to predictions of diffractive dijet photoproduction to 
the corresponding ratio in DDIS [18]. In this way several experimental systematic uncertainties 
cancel and theoretical uncertainties can be reduced. 
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In the present paper new measurements of diffractive dijet cross sections in DIS and photo¬ 
production are presented. The data were collected in the years 2006 and 2007 with a total in¬ 
tegrated luminosity of 30 pb 1 for diffractive photoproduction and 50 pb 1 for diffractive DIS. 
For the identification of diffractive events a proton detected in the VFPS is required. The results 
are compared to NLO QCD calculations. 


2 Kinematics 


Figures |T|( a) and (b) show leading order diagrams of direct and resolved diffractive dijet pro¬ 
duction in ep interactions. The relative contributions of these two components depend on the 
virtuality of the exchanged photon such that at high virtualities the direct process is dominating. 
The incoming (scattered) positron four-momentum is denoted as k (A/), the four-momentum of 
the virtual photon emitted from the positron as q = k — k'. The four-momentum of the incoming 
(outgoing) proton is P ( P '). The kinematics of the ep scattering process can be described by 

s={k + P) 2 , Q 2 = —q 2 , 

where s is the square of the centre-of-mass energy of the collision, Q 2 is the photon virtuality 
and y the inelasticity of the process. With P x being the four-momentum of the hadronic final 
state excluding the leading proton (see figure [T]), the inclusive diffractive kinematics is described 
by the additional variables 

M x — Px, x jp = q '^~ p P, \ t={P-P')\ (2) 

where M x is the invariant mass of system X, xu> corresponds to the longitudinal momentum 
fraction lost by the incoming proton and t is the four-momentum transfer squared at the proton 
vertex. 


For diffractive dijet production additional invariants are introduced. With denoting the four- 
momenta entering the hard sub-process from the photon and from the pomeron side as u and 
v, the longitudinal fractions of the photon and of the pomeron momentum entering the hard 
sub-process, x 7 and z F , are defined as 


x 7 = 


P ■ u 
P ■ q 


and 


Z]p — 


q ■ v 


q-{P-P'Y 


(3) 


respectively. 


In leading order, the invariant mass of the dijet system M 12 is equal to the centre-of-mass 
energy of the hard sub-process 


M\2 = (u + Vf. 


(4) 
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Figure 1: Leading order diagrams of the direct a) and resolved b) diffractive dijet production. 


3 Factorisation in Diffractive Dijet Production 


In the QCD factorisation approach the diffractive dijet cross section is given by the convolution 
of partonic cross sections da with diffractive parton distributions f[j p : 


d a(ep —> e + 2 jets + X' + p) — J dt J dxjp J d zjp 

d(J e i—y2 jets(®j Id' ) X fi/pizFiVFiXPit)- 


(5) 


Here, the hadronic system X' corresponds to what remains of the system X after removing the 
two jets. The integrals extend over the accepted phase space. The sum runs over all partons i 
contributing to the cross section, s ~ xjpzpys — Q 2 is the sub-process invariant energy squared 
and up and //\> denote the factorisation and renormalisation scales, respectively. 

In the photoproduction region the exchanged photon may dissociate into a low mass non- 
perturbative hadronic system due to its low virtuality (figure [TJ d) and a photon parton distribution 
function ( 7 PDF) is introduced. The cross section for this resolved photon process is given by 


dcr(ep —> e + 2 jets + X' + p) = j dt J dxjp j 




ip 


d V f"//e{y) I dx 7 f j/l {x 1 ^ Z F ) X dd-^2jets(s,^R,/^F) X A(F) X]p, f), 


( 6 ) 


where / 7 / e is the Weizsacker-Williams equivalent photon flux [23, [24] integrated over the mea¬ 
sured Q 2 range and /j / 7 are the parton distribution functions in the photon ( 7 -PDF). In this 
case, the centre-of-mass energy of the hard subprocess is approximated by s ~ x^xpzn-ys. As 
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default, the GRV [25Q 7 -PDFs are used to describe the structure of resolved photons. The AFG 
7 -PDF set is also studied. 


For the diffractive proton parton densities, the HI2006 Fit-B DPDF set [4] is used. This 
parametrisation was obtained from a QCD fit in NLO accuracy to inclusive DDIS data. In this 
fit a proton vertex factorisation [27] is assumed in which the xp and t dependencies of the 
DPDFs factorise from the dependencies on (i F and zp such that 

fi/ p ( z ip-:E fi x ipi t) — hp/p{ x iPi t) fi/p{zp, p F ) + rip fp/ p (x jp, t) fi/p{zp, p F ). (7) 

The pomeron flux factor fp/ p {xp,t) was parametrised in [0] as suggested by Regge models 
[28]. For xp 0.01 a small additional contribution from sub-leading reggeon ( R) exchange 
described by the second term in ([7]) was taken into account, where rip ~ 10 -3 is the normalisa¬ 
tion factor of the reggeon contribution [ 0 ] • 


4 NLO QCD Calculations 


Theoretical calculations of dijet production in next-to-leading order were performed in the 7 p 
regime using the the FKS program [29] and in DIS using NLOJET++ [30, 31]. Both programs 
were adapted lfl8ll for hard diffraction. The NLO calculations for photoproduction are consistent 


with calculations performed by Klasen and Kramer [32-34]. Similarly, the DDIS predictions 
were checked using the independent package DISENT NLO [35]. 


The NLO calculations are performed with the number of flavours fixed to 5 and the QCD 
scale parameter set to A 5 = 0.228 GeV, corresponding to a 2-loop as{M z ) of 0.118. The 
renormalisation and factorisation scales are set to be equal and are calculated from the average 
jet transverse energy (E^ et ) = (E^ etl + E^ et2 )/ 2 and the momentum transfer Q 2 as ji F = 
li F = (E^ ct ) 2 + Q 2 . Lor photoproduction, Q 2 is set to zero. The sensitivity of the NLO 
predictions to the scale choice is studied by varying the scale up and down by a factor of two. 
An alternative definition of the scale n 2 R — n 2 F — (E^ et1 ) 2 + Q 2 / 4, based on the leading jet 
transverse energy E^ etl , is also studied. 


5 Experimental Procedure 


5.1 The HI Detector 


A detailed description of the HI detector can be found elsewhere [36-38]. Here only the de¬ 
tector components most relevant to the present analysis are briefly described. A right-handed 
coordinate system is employed with the origin at the nominal ep interaction point and with the 
positive z-axis pointing in the proton beam direction. The x-axis is pointing along the horizon¬ 
tal direction to the centre of the HERA ring. The pseudorapidity 7 = — In tan | is calculated 
using the polar angle 9 measured with respect to the proton beam direction. 

The interaction point is surrounded by the central tracking detector (CTD), which consists of 
a set of concentric drift chambers supplemented by silicon detectors [39] located inside the drift 
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chambers. Charged particle trajectories are bent by a 1.15 T homogeneous solenoidal magnetic 
field. The region in pseudorapidity covered by the CTD is —2.0 < r] < 2.0 and the trans¬ 
verse momentum resolution is a {pt)/pt — 0.002 p T /GeV © 0.015. A multi-wire proportional 
chamber at inner radii (CIP) is mainly used for triggering [40]. The forward tracking detector 
supplements the CTD track reconstructions in the region 7° < 6 < 25°. 

Scattered positrons in the rapidity range —4 < r/ < —1.4 are measured in a lead / scintillating- 
fibre calorimeter, the SpaCal [38], with energy resolution 7%/ \JE/GcV © 1%. 


The central and forward tracking detectors are surrounded by a finely segmented Liquid 


Argon (LAr) calorimeter [41] situated inside the solenoidal magnet and covering the pseudo¬ 
rapidity region —1.5 < rj < 3.4. Its resolution was measured in test beams 1142114311 and is 
11%/ \JE /GcV © 1% and 50%/^/77/GeV © 2% for electromagnetic and hadronic showers, 
respectively. The hadronic energy scale is known within 2% for this analysis a. 

The ep luminosity is determined online by measuring the event rate of the Bethe-Heitler 
bremsstrahlung process, ep —> epp, where the photon is detected in a calorimeter located close 
to the beam pipe at z = —103 m [36]. The overall integrated luminosity normalisation is 
determined using a precision measurement of the QED Compton process [45J]. 


5.2 Very Forward Proton Spectrometer 

The Very Forward Proton Spectrometer (VFPS) consists of two Roman pots located 218 m and 
222 m from the interaction point in the forward direction. It allows for a measurement of protons 
with energies between 895 and 912 GeV (0.008 < Xp < 0.028) and with transverse momenta 
up to about 0.8 GeV (\t\ < 0.6 GeV 2 ) (J. 

The VFPS complements the Forward Proton Spectrometer (FPS) HI]. The FPS has a wider 
acceptance in scattered proton energy (xp < 0.1) but has only limited geometrical acceptance 
in the azimuthal angle of the scattered proton (figure [2]). In particular at small \t\ < 0.2 GeV 2 , 
the VFPS acceptance is much better than for the FPS. More than 70% of the diffractive events 
have \t | smaller than 0.2 GeV 2 . 

The Roman pots, which are moved close to the beam as soon as the beam conditions are 
sufficiently stable, are equipped with detectors made of several layers of scintillating fibers with 
photomultiplier readout. The sensitive detector areas are covered by scintillator tiles, the signals 
of which are used as a trigger. The VFPS has high track efficiency (~ 96%) and low background 
contamination (~ 1%). 


5.3 Kinematic Reconstruction 


The observable xp is reconstructed by the VFPS from the relative distance and angle between 
the track reconstructed between the two stations and the beam and can be expressed as 


Xp 



( 8 ) 
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Figure 2: Beam envelope [0] as a function of the distance s to the HI vertex in the x projection, 
for the p beam and diffractive protons at xjp = 0.02 and \t\ = 0.01, 0.1 and 0.5 GeV 2 . The 
locations of FPS/VFPS stations are indicated by the vertical lines. 


where E' is the energy of the leading proton in the VFPS and E p is the proton beam energy. 

The quality of the reconstruction of ./■)/!' 1 ’ s was checked using an event sample of elasti¬ 
cally produced p mesons, ep —> epp. The x F variable reconstructed from the p decay tracks 
detected in the CTD is compared to x]p FPS determined by the VFPS stations. The resulting 
_ x Pp distributions are found to be in good agreement with Monte Carlo simulations [46]. 
The resolution of x)p PS was determined to be equal to 0.0022 [0]. 


The invariant mass M x of the system X is calculated from all hadronic objects in the main 
HI detector: 



(9) 


The hadronic final state (HFS) is reconstructed using an energy flow algorithm which combines 
information from the trackers and calorimeters by avoiding double-counting of energies [47, 


48]. 


Jets are reconstructed from the hadronic final state objects using the longitudinally invariant 
kx -jet al gor ithm [49Q with a jet distance parameter R = 1.0 as implemented in the FastJet 
package [50]. The massless p^-recombination scheme is used. The jet finding algorithm is 
applied in the 7 *p frame. In photoproduction this frame is identical to the laboratory frame up 
to a Lorentz boost along the beam axis. 
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The jet properties are studied in terms of the transverse energy of the leading jet E^ etl in 
the 7 *p frame, of the invariant mass of the dijet system M 12 and of the pseudorapidity variables 
| Ar/i ets | and (// Jets ) defined in the laboratory frame, where 


Ml = 

(J (1 > +J< 2) ) 2 , 

(10) 

|A/ 7 jets | = 

7*“ - rf a2 \, 

(ID 

W ets ) = 

\ 7 etl + v et2 ) ■ 

(12) 


In these definitions, J 0) and denote the four-momenta of the two leading jets. 


5.3.1 DIS 


For DIS events the polar angle 9' e and energy E' e of the scattered positron are measured in the 
SpaCal calorimeter. The kinematic reconstruction method introduced in [511] is used 


V = VDA+ vl - vIa, Q 2 = (13) 

tan 2 ^ 

This method interpolates between y e determined from the scattered positron alone at larger 
inelasticity and ijda determined using the double angle method at low y. 

The variable zp is calculated as 


„obs 


Q 2 + m 2 2 
Q 2 + Ml ' 


(14) 


5.3.2 Photoproduction 


In the 7 p regime the scattered positron leaves the interaction undetected. Therefore the inelas¬ 
ticity y is reconstructed from the hadronic final state 


V = 


Pz,i) 

2 E e 


where E e is the initial positron beam energy. 

The observables xp and zp are calculated from the hadronic final state X as 


(15) 


x 


obs 

7 


S^zGjets (Ei 
^2iex(Ei 



and zp 


Szgjets(-^ + E Zt j) 


(16) 


where the sums in the numerators run over the leading and the sub-leading jet, whereas the sums 
in the denominators include all objects of the reconstructed hadronic final state. 
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5.4 Event Selection 


The analysis is based on a data sample corresponding to an integrated luminosity of 30 pb 1 for 
photoproduction and 50 pb -1 for DIS collected with the HI detector in the years 2006 and 2007 
with proton and positron beam energies of 920 GeV and 27.6 GeV, respectively. The events are 
triggered on the basis of a coincidence of VFPS signals from both stations, together with condi¬ 
tions on the charged track transverse momenta and track topology in the HI main detector [52]. 
The trigger efficiency, calculated using events collected with an independent trigger condition, 
was found to be about 80% with negligible dependence on kinematic quantities. This efficiency 
is well reproduced by the HI trigger simulation after correcting for an overall normalisation 
difference of 5%. For the DIS analysis the integrated luminosity is increased using the fact that 
for most of the DIS events also another trigger based on signals in the SpaCal has fired. Only 
events with a VFPS track in a fiducial volume of high efficiency are selected [Q]. The recon¬ 
structed ^-coordinate of the event vertex is required to be within 30 cm of the mean ^-position 
of the interaction point. 


The random overlap of ep events with beam-halo protons detected in the VFPS can consti¬ 
tute a possible background to the VFPS diffractive data sample. In such background events the 
detected proton typically has a small energy loss, not compatible with the energy loss expected 
from the energy deposited in the main HI detector. The relative energy loss of the proton de¬ 
tected in VFPS, x^ PS , is thus required to be at least 60% of x p ) measured in the HI dctcctoiQ, 
x]p PS /x]p > 0.6. In addition, .r 1 /, 1 is required to be smaller than 0.04. The remaining back¬ 
ground contamination after applying the above cuts is estimated from data by overlaying events 
without VFPS activity with VFPS signals recorded independently of any detector activity and 
is found to be less than 1 % [ 2 ]. 


The scattered positron candidate of an event is identified as the electromagnetic cluster with 
the highest transverse momentum being well isolated and having a minimum energy of 8 GeV. 
If such a candidate is absent the event is defined as photoproduction. 


For the selection of DIS events in this analysis the positron candidate is required to be 
detected in the SpaCal. The energy E' e and polar angle 9' e of the scattered positron are deter¬ 
mined from the SpaCal cluster and the interaction vertex reconstructed in the CTD. In order 
to improve the background rejection, additional requirements on the transverse cluster radius 
and lower limit to the positron energy are imposed 15311 . The quantity J2i(Ei ~ P z ,i) summed 
over all HFS particles and the scattered positron, is required to be in the range 35-75 GeV. For 
fully reconstructed neutral current DIS events this quantity is expected to be twice the positron 
beam energy (55.2 GeV) but is expected to be lower for photoproduction background where the 
scattered positron escapes undetected. Radiative events where a photon is emitted along the 
direction of the incident positron beam, also have a reduced ~ -Pz.i)- 


The leading and the sub-leading jets are required to have transverse energies E, 


*jetl 
T 


> 


5.5 GeV and E * |e12 >4.0 GeV, respectively. These cuts are asymmetric in the transverse en¬ 
ergy to restrict the phase space to a region where NLO QCD calculations are reliable [54, 55]. 
An event is rejected if one of these two jets is outside of —1 < rf et1 ' 2 < 2.5. Events with z°p s 


'The variable xf, 1 is calculated as 
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Photoproduction DIS 

Event kinematics 

Q 2 < 2 GeV 2 4 GeV 2 < Q 2 < 80 GeV 2 

0.2 < y < 0.7 

Diffractive phase space 

0.010 < x P < 0.024 
\t\ < 0.6 GeV 2 

Zjp < 0.8 

Jet phase space 

£; jetl > 5.5 GeV 

E^ et2 >4.0 GeV 
—1 < r/J et1 ’ 2 < 2.5 


Table 1: Phase space of the diffractive dijet VFPS measurement for photoproduction and deep- 
inelastic scattering. 

above 0.8 are excluded to improve the reliability of the comparison between data and theoretical 
predictions, since the DPDF are determined with a similar zjp restriction. 

The DIS events are selected with photon virtualities 4 GeV 2 < Q 2 < 80 GeV 2 . Both data 
samples are restricted to a common y range 0.2 < y < 0.7. In table Q] the definitions of both 
analysis phase spaces are summarised. The photoproduction and DIS data samples contain 
3768 and 550 events, respectively. In addition to the event selection summarised in table |Tj 
an event selection is performed extending the phase space in all kinematic variables and other 
selection requirements to obtain events for an adequate description of migrations at the phase 
space boundaries. 


5.5 Monte Carlo Simulations 

5.5.1 Correction to the Data 


The Monte Carlo (MC) simulation method is used to correct the data for effects of detector 
acceptance, resolution and detector inefficiencies. All MC samples are passed through a de¬ 
tailed HI detector simulation based on the GEANT program [56] and are subjected to the same 
analysis chain as is used for the data. 


Diffractive dijet photoproduction and DDIS events were generated using the RAPGAP MC 
generator [57]. This generator is based on leading order (LO) parton level QCD matrix elements 
with a minimum transverse momentum of the outgoing partons of = 1.7 GeV. Higher or¬ 
ders are mimicked by initial and final state leading logarithm parton showers. Fragmentation 
is accounted for using Lund string model [58] as implemented in Pythia MC generator [59jJ. 
The HI2006 Fit-B DPDF set [4] is used in RAPGAP to describe the density of partons in the 
diffractively scattered proton. In photoproduction a resolved photon contribution is simulated 
using the GRV-LO photon distribution function |f25l] . In addition to a pomeron exchange con¬ 
tribution also a sub-leading reggeon contribution is included, corresponding to about ~ 2% of 
the total cross section. In order to describe the data sufficiently well reweighting functions are 
applied in z$ s , xn> and t. The reweighting is different for zp and DIS. 
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5.5.2 Correction to Theoretical Models 


For comparison of QCD calculations with the diffractive measurements, it is necessary to con¬ 
vert the calculated NLO parton level cross sections to the level of stable hadrons by evaluating 
effects due to hadronisation, fragmentation and the influence of pomeron or photon remnants. 
The RAPGAP MC generator is used to compute the required hadronisation correction factors 
for the diffractive dijet calculations. These factors are defined for each measured data point by 


1 + <5* — 


a; 


hadr 


(7 


part ' 


(17) 


where the o diadr (of ait ) are the bin-integrated MC cross sections at hadron level (parton level) in 
a given bin i. They reduce the predicted NLO parton level cross sections by typically ~ 9% in 
photoproduction and enhance the cross sections by typically ~ 2% in DIS. In photoproduction 


the hadronisation correction factor is particularly large at the second highest x° bs bin, where 
contributions with x° bs ~ 1 at parton level migrate to lower values due to hadronisation effects. 
The hadronisation corrections have uncertainties of 3% [17]. The hadronisation corrections 
determined here are applicable to NLO QCD predictions, since a good agreement in shape of 
the parton level predictions of the MC to the NLO calculations is observed. 

In the DIS analysis, the RAPGAP MC generator is also used to correct the measured data 
for QED radiation effects. The radiative corrections are defined as 


1 A* — 
- 1 - 1 


c7, 


nrad 


cr; 


rad ’ 


(18) 


where cr[ ad (<rf rad ) denote the bin integrated cross sections obtained from RAPGAP when run 
with (without) simulating QED radiation. The term <5* ad is on average compatible with zero with 
a standard deviation of 4% within the phase space. Radiative corrections in photoproduction 
are found to be negligible. 


For the comparison with the measurement, the NLO QCD predictions are scaled down by a 
factor of 0.83 [60] to account for the contributions from proton dissociation (My < 1.6 GeV) 
absent in the current analysis but included in the extraction of the HI2006 Fit-B DPDF set from 
the inclusive data [0]. 


5.6 Cross Section Measurement 


In order to correct for detector effects, the dijet cross sections are calculated at the level of 
stable hadrons using a matrix unfolding method En, [62]. The detector response is described by 
a matrix A determined from the RAPGAP simulation. It relates the expected vector of event 
counts, (i/rec), to the true event count vector, x tTue , on the level of stable hadrons via the formula 

(i/rec) AXt rue . 


In order to control migrations at the phase space boundaries also the neighbouring parts of 
the analysis phase space are taken into account. Of these, the migrations caused by events in 
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which jets have low transverse momenta E T , high xjp or low y are most important. Similar 
unfolding techniques have been applied in other jet-based analyses Iild.llll. l631. 

An estimator of the true-level event count x trU e is obtained by minimising a y 2 function (fl9l) 
with respect to x tIue 


1 


X^ XA~b^~ Xl r> (Z/rec -^-^true) V (?/rec A't'true)(■t'true -^b) L L^true •t'b)) (19) 


T t T i 


where the matrix V is the covariance of data. The term y^ is a measure of the agreement 
between A.x true and y Tec , where y rec is the vector of events counts after subtracting residual 
background contributions. The regularisation term r 2 y| suppresses large fluctuations of :H tnie . 
The type of the regularisation is defined by the matrix L. In this paper, L is set to the unity 
matrix. The vector defines a bias for the regularisation term, taken from the RAPGAP 
prediction. The value of the regularisation parameter r is chosen using the L-curve method as 
described elsewhere 10 . 


The bin-integrated cross section for each data point is given by 


« = % (1 + C) ( 20 ) 

where £ is the integrated luminosity of the data. The radiative corrections <5 rad are non-zero only 
for the DIS case. 


5.7 Systematic Uncertainties of the Measured Cross Section 

For each source of systematic uncertainty, a separate response matrix A is filled and the dif¬ 
ference to the nominal matrix A is propagated through the unfolding procedure. All these 
individual contributions of systematic uncertainties are then added in quadrature for each bin to 
obtain the total systematic uncertainty. The following systematic effects are studied: 

VFPS calibration The primary source of the VFPS systematic uncertainties is related to an 
uncertainty of the x and y global track coordinate^ with respect to the beam. The actual 
beam position is measured with help of a beam position monitor [2|, 64] which has a 
precision of 160 yn i in x and 120 //m in y. The horizontal coordinate x has an additional 
uncertainty originating from the VFPS calibration procedure, tied to the reconstruction of 
xjp in the main HI detector. The resulting /(-coordinate uncertainty is 250 //m. 

The time variation of the beam-tilt in x and y introduces an uncertainty of 8 //rad for the 
x-tilt and 6 //rad for the y-tilt. 

More details on the VFPS reconstruction and its precision are given in [2i]. In total, all 
sources of the VFPS uncertainties affect the integrated cross section by 5.5% in yp and 
typically 3.7% in DIS. 

2 The global track coordinates are reconstructed by linking the local tracks of the two VFPS stations. 
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Positron reconstruction In the DIS analysis the uncertainties of the measured positron energy 
E' e (1%) and angle 9' e (1 mrad) in the SpaCal calorimeter lead to an uncertainty of the total 
cross section of 0.4% and 0.7%, respectively. 


Energy scale The uncertainty of the hadronic final state energy calibration is 2% [44]. It affects 
the total cross section by ±7.6% for photoproduction and by ±6.1% for DIS. 


Model uncertainties The influence of the MC model used to unfold the cross sections is stud¬ 
ied by varying the kinematic distributions of the RAPGAP MC generator within cer¬ 
tain limits while maintaining an acceptable description of the data. For this purpose 
the shape of the kinematic distributions in E^ etl , xjp, z n >, x 7 , y, t and Q 2 are altered 
by applying multiplicative weights of (E^ etl ) ±0A , x|> 0 ' 2 , z% 03 , x± 0A , y ±0 3 , e ±f and 
(i Q 2 ±0.1 GeV 2 ) ±0 - 2 , respectively. The largest resulting uncertainties arise from vari¬ 
ation of the shape in t (4.5% in yp and 3.3% in DIS) and E^ ctl (3.5% in yp and 3.0% in 
DIS). The integrated cross section uncertainty due to model dependence is 7% in yp and 
5% in DIS. 


Normalisation uncertainties The following sources of systematic normalisation errors are 
considered: 


• The VFPS track reconstruction efficiency is known to within 2.5% 


10 . 


• The VFPS background originating from interactions of beam particles with the 
residual gas, producing a proton signal in the VFPS in accidental coincidence with a 
dijet event in the main HI detector is less than 1% and is treated as a normalisation 
uncertainty 0. 


The integrated luminosity of the VFPS triggered data is known to within 3% [450 - 


• The trigger efficiency has an uncertainty of 5%. 


The resulting total normalisation uncertainty amounts to 6%. 


Figure [3] displays the distributions of the x- and //-coordinates of global tracks in the VFPS, 
xjp as measured by the VFPS, Q 2 for the DIS selection and the jet variables E^ etl and (% ets ) in 
comparison to the MC distributions after reweighting and normalising to the data. In all cases 
the data are well described in shape within systematic errors. 


6 Results 

6.1 Integrated Photoproduction and DIS Cross Sections 

The integrated e + p diffractive dijet cross sections in the yp and in the DIS regime measured 
in the kinematic range defined in table[T| are presented together with NLO QCD and RAPGAP 
predictions in table [2] The total theoretical uncertainty is calculated by using the sign improved 
quadratic sum of DPDF eigenvectors a, scale and hadronisation uncertainties. In the DIS 
regime, the theoretical expectation agrees with the measurement within uncertainties. This 
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PHP 

DIS 

Data [pb] 

237 ±14 (stat) ±31 (syst) 

30.5 ±1.6 (stat) ±2.8 (syst) 

NLO QCD [pb] 

430 tH 2 ( scale ) tei (DPDF) ±13 (hadr) 

28.3 tl l 4 4 (scale) tit (DPDF) ±0.8 (hadr) 

RAPGAP [pb] 

180 

18.0 

Data/NLO 

0.551±0.078 (data)l^ (theory) 

1.08±0.11 (datajlg;^ (theory) 


Table 2: Integrated e + p diffractive dijet cross sections in 7 p and DIS compared to NLO QCD 
calculations using the HI2006 Fit-B DPDF set. The measured cross sections are gresented with 
statistical and systematical uncertainties. For the theoretical predictions the uncertainties from 
scale variations, from the H12006 Fit-B DPDFs and from the hadronisation corretions are given. 
The predictions from RAPGAP are also shown. The ratios data/NLO are given in the last row. 


confirms the observation made in previous measurements [[ 8 - IT]. In contrast, the integrated 
diffractive dijet cross section in photoproduction is overestimated by the NLO QCD theory by 
almost a factor of two, with considerable theory uncertainty. This observation agrees with the 
results of previous HI analyses in a similar ki nematic range lfl7l Hill. based on different data 
sets and using different experimental techniques to select diffractive events. To conclude, the 
integrated NLO QCD cross section predictions are in disagreement with three independent HI 
measurements of diffractive dijet photoproduction. The MC RAPGAP, based on leading order 
matrix elements and parton showers, fails to describe the integrated cross sections both in DIS 
and in photoproduction. 


6.2 Diffractive Dijet Production in DIS 


The measured differential DIS cross sections as a function of Zp s , xjp, y, Q 1 are given in table [4] 
and are shown in figure |4] together with the NLO QCD predictions. In table [5] and figure [5] the 
differential cross sections in DIS are shown as a function of E^ ctl , (// Jels ), Ar/ Jets and M x ■ The 
NLO QCD predictions are in good agreement with the measurements within data and theory 
uncertainties. 


The shapes of the NLO predictions are tested using the ratio of data to prediction. A some¬ 
what different shape is observed for data and theory as a function of Q 2 , however the deviations 
are covered by the uncertainties. Resolved photon [65] and higher twist contributions [661 are 
expected to change the DIS cross sections at small Q 2 . The predicted shape in y also differs from 
the observation, such that at high y smaller cross sections are predicted than observed. Similar 
shape deviations in Q 2 and y have also been observed in a recent measurement of diffractive 
dijet production based on a large rapidity gap selection |Q. The cross section as a function 
of E^ ctl is observed to be slightly harder than predicted by theory, although still in agreement 
within uncertainties. 


6.3 Diffractive Dijet Production in Photoproduction 

The measured differential cross sections in the 77 ;-regime are given in table [ 6 ] and shown in 
figure [ 6 ]as a function of z^ s , xjp, y, x° bs together with the NLO QCD calculations. The differ- 
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ential cross sections for the variables E^ etl , (r/ je1s ), ] Ay iRts j and M x are given in table |7] and 
shown in figure [7] The relative statistical uncertainties in photoproduction are in most cases 
smaller than in the case of deep-inelastic scattering. The NLO QCD predictions agree well with 
the measured distributions in shape but overestimate the dijet cross sections in normalisation, 
as already discussed for the integrated cross sections. In particular there is no significant depen¬ 
dence of the data to theory ratio on the variables z'jjl *, x° bs and E^ ctl which are sensitive to the 
DPDF and to the presence of a diffractive exchange remnant. These results are in qualitative 
agreement with previous HI measurements [ 17l 18]. Using the AFG [260 photon PDF as an 
alternative the predicted integrated cross section is reduced by 6%. As visible in figures |6] and 
0 the shapes of the distributions depend only little on the choice of the photon PDF. 


6.4 Comparison of Dijet Cross Sections in Diffractive Photoproduction 
and DIS 

The conclusions made in previous sections about the normalisation problems of the NLO calcu¬ 
lations in diffractive photoproduction suffer from large theoretical uncertainties. This situation 
is summarised in figure [8] where the ratio of observed cross section to expectation is shown 
as a function of Q 2 , also including an extra bin for the cross section in the photoproduction 
regime, 0 < Q 2 < 2 GeV 2 . No significant deviation from unity is observed for the suppression 
fact oi0 as a function of Q 2 in the DIS regime, whereas the NLO calculation fails to describe 
the measurement in the photoproduction region. For comparison, also the ratio of the RAPGAP 
prediction to the NLO calculation is shown. RAPGAP is off in normalisation and predicts a 
shape in Q 2 which differs from the NLO calculation. 

In a refined method for studying deviations of the NLO QCD predictions from photopro¬ 
duction data the cross sections measured in the yp-regime are divided by the corresponding 
cross sections in DIS. In such ratios most of the data systematic uncertainties are reduced, with 
the exception of the model uncertainties which are uncorrelated between yp and DIS. Similarly, 
theoretical uncertainties cancel to a large extent. This is true for the DPDF uncertainties as well 
as for scale variations, if the NLO QCD scales are varied simultaneously for photoproduction 
and DIS. The hadronisation corrections, however, are taken to be uncorrelated between DIS 
and photoproduction, such that they amount to about y/2 x 3% in the ratio of the integrated 
cross section. The resulting cross section ratios of photoproduction to DIS are summarised in 
table [3] The double-ratio of photoproduction to DIS, data to NLO, is also given and shown 
in figure 0 Due to the reduced theoretical uncertainty the double ratio deviates significantly 
from unity indicating that factorisation does not hold in diffractive dijet photoproduction with 
respect to the same process in DIS. This statement is valid within the theoretical framework 
applied in this paper and under the assumption that the scale must be varied simultaneously for 
the DIS and yp calculations, which leads to cancellations of the respective uncertainties in the 
ratio. Higher order corrections may change this picture. As an estimate of the possible size 
of such corrections the difference between leading-order and NLO calculations scaled by ay/2 
may be taken, which amounts to 5%. When changing the photon PDF from the GRV PDF set 
to the AFG PDF set a rise in the double ratio of 6% is observed. Using p 2 = E^ etl2 + Q 2 /4 

3 The suppression factor is defined as a ratio of data and NLO QCD cross section. 
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Ratio of photoproduction to DIS 

Data 

NFO QCD 

7.78 ±0.60 (stat) ±1.14 (syst) 

15.21 1°;°° (scale) 1°; 22 (DPDF) ±0.65 (hadr) 
14.22 with AFG 7 PDF 

14.17 with scale fj, 2 = (£’y ? et1 ) 2 + Q 2 /4 

Data/NFO 

0.511 ±0.085 (data) +°;° 22 (theory) 

0.547 with AFG 7 PDF 

0.548 with scale /x 2 = (E^ et1 ) 2 + Q 2 /A 


Table 3: Ratio of integrated e + p diffractive dijet cross sections for () 2 < 2 GeV 2 (photoproduc¬ 
tion) to Q 2 > 4 GeV 2 (DIS). Listed are the ratios for data and for the NLO calculation including 
two variants. The data and NLO uncertainties are indicated. The double-ratio of data to NLO 
and its uncertainties are also given. 


as the scale choice leads to an increase of the double ratio by 7%. The observed suppression 
agrees with previous HI results EM n is worth mentioning that the suppression is now 
measured at HERA both in processes with an identified leading proton and in processes with a 
large rapidity gap selection, so possible contributions from proton-dissociative processes alone 
are excluded as an explanation. 

Possible shape dependencies of the suppression are studied using cross section ratios of pho¬ 
toproduction to DIS differential in the variables |A 77 jets |, y, zp and Ej! et 1 , as given in table [ 8 ] 
The data ratios as a function of | A 7 / jets | and y are shown in figure ITOl together with predictions 
from NLO QCD and RAPGAP. The measured shapes are not described well, but the limited ex¬ 
perimental precision does not allow for strong conclusions to be made. The ratios as a function 
of zp and E^ ctl are shown in figure Qj] Within uncertainties the corresponding double ratios 
are constant throughout the measured z F and E^ ctl ranges. 


7 Summary 


Diffractive dijet production is measured in photoproduction and deep-inelastic scattering in the 
same kinematic range 0.2 < y < 0.7 and 0.010 < xp < 0.024 for jets with E^ etl > 5.5 GeV, 
E^ et2 > 4.0 GeV and with limits on the photon virtuality Q 2 < 2 GeV 2 for photoproduction and 
4 GeV 2 < Q 2 < 80 GeV 2 for DIS. For the leading proton detection, the HI Very Forward Proton 
Spectrometer is used for the first time, such that the diffractive sample is free of background 
from low-mass proton dissociative states. 

In DIS, diffractive dijet production is well described within the experimental and theoretical 
uncertainties by the NFO calculations based on the HI2006 Fit-B diffractive parton densities 
of the proton. Within uncertainties, the QCD factorisation assumptions made for the NFO 
calculation are confirmed in this process. This result is consistent with previous HI and ZEUS 
measurements and the new data may be used in future DPDF fits. 
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In photoproduction, next-to-leading order (NLO) calculations based on the HI2006 Fit- 
15 diffractive parton densities overestimate the measured total cross sections, thus confirming 
previous HI measurements, where the Large Rapidity Gap method for the identification of 
diffractive events was used. The shapes of the differential cross sections are described within 
the experimental and theoretical uncertainties. There is no hint of dependence of the observed 
suppression on the variable x° bs . 

In order to profit from cancellations of theoretical uncertainties, ratios of photoproduction to 
DIS cross sections and double ratios of data to NLO are analysed. Integrated over the analysis 
phase space the double ratio is found to be 0.51 ± 0.09. Following this, within the theoreti¬ 
cal framework based on diffractive parton densities, factorisation is broken in diffractive dijet 
photoproduction. This observation is in agreement with previous HI measurements, where 
complementary experimental methods have been used. Contributions from proton dissociative 
processes present in the previous analyses but absent here are ruled out as a cause of the ob¬ 
served suppression. The differential measurements of cross sections and cross section ratios 
in DIS and photoproduction provide stringent tests of the theory both in normalisation and in 
shape. 
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Table 4: Integrated cross section and bin averaged hadron level differential cross sections as a function of the variables zp, xp, y and Q 2 for 
diffractive dijet DIS in the phase space detailed in tabled] For each data point, the statistical (<5 stat ) and systematic (d sysl ) uncertainties and the 
corresponding correlation coefficients (p stat , p sys t ) are given. The hadronisation correction factors (1 + <5 ha dr) applied to the NLO calculations 
and the radiative corrections (1 + 5 ra d) are also listed. The overall normalisation uncertainty of 6% is not included in the table. 
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Table 5: Bin averaged hadron level differential cross sections for diffractive dijet DIS as a function of the variables E^ et , Mx, | Ar/ Jels and 
(r/' ets ) in the phase space detailed in tabled] For each data point, the statistical (<5 stat ) and systematic (5 syst ) uncertainties and the corresponding 
correlation coefficients (p stat , p syst ) are given. The hadronisation correction factors (1 + r) hadr ) applied to the NLO calculations and the the 
radiative corrections (1 + <5 rac j) are also listed. The overall normalisation uncertainty of 6% is not included in the table. 
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Table 6: Integrated diffractive dijet ep cross section and bin averaged hadron level differential diffractive dijet ep cross sections as a function 
of the variables z P , xjp, y and x 1 for the dijet photoproduction kinematic range in the phase space detailed in tabled] For each data point, 
the statistical (<5 stat ) and systematic (5 syst ) uncertainties, the corresponding correlation coefficients (p stat , p syst ) and the hadronisation correction 
factors (1 + <5 hadr ) applied to the NLO calculations are given. The overall normalisation uncertainty of 6% is not included in the table. 
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Table 7: Bin averaged hadron level differential diffractive dijet ep cross sections in the photoproduction kinematic range as a function of the 
variables E^ et , M x , A;/' cts and (// |ets ) in the phase space detailed in table Q] For each data point, the statistical (5 stat ) and systematic (<5 syst ) 
uncertainties, the corresponding coefficients (p stat , p syst ) and the hadronisation correction factors (1 + c) hadr ) applied to the NLO calculations 
are given. The overall normalisation uncertainty of 6% is not included in the table. 
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Table 8: Ratios of differential diffractive dijet ep cross sections, measured in photoproduction, to measurements in DIS as a function of 
the variables |A// ets |, y, zjp and E^ tU in the phase space detailed in table [H For each data point, the statistical (5 stat ) and systematic (5 syst ) 
uncertainties, the corresponding correlation coefficients (p stat , p sys t) and the hadronisation correction factors (1 + <5 hadr ) applied to the NLO 
calculations are given. 
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Figure 3: Comparison of the photoproduction data (dots) and DIS data (triangles) with the 
reweighted RAPGAP MC simulation (solid line) as a function of coordinates x and y in VFPS, 
xip, Q 2 , Ejt ei] and (r/ je1s ). The systematic uncertainties are shown as bands on the histograms. 
For better visibility the (// Jet ' s ) DIS distribution is multiplied by a factor 3. 
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Figure 4: Diffractive dijet DIS cross sections differential in zp, xp, y and Q 2 . The inner error 
bars represent the statistical errors. The outer error bars indicate the statistical and systematic 
errors added in quadrature. The overall normalisation uncertainty of 6% is not shown. NLO 
QCD predictions based on the HI2006 Fit-B DPDF set, corrected to the level of stable hadrons, 
are shown as a white line. They are scaled by a factor 0.83 to account for contributions from 
proton-dissociation which are present in the DPDF lit but not in the data. The inner, light 
shaded band indicates the size of the DPDF uncertainties and hadronisation corrections added 
in quadrature. The outer, dark shaded band indicates the total NLO uncertainty, also including 
scale variations by a factor of 0.5 to 2. For each variable, the cross section is shown in the upper 
panel, whereas the ratio to the NLO prediction is shown in the lower panel. 
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Figure 5: Diffractive dijet DIS cross sections differential in E^ ctl , M x , A;/ Jets and (r/ jets ) . The 
inner error bars represent the statistical errors. The outer error bars indicate the statistical and 
systematic errors added in quadrature. Further details are given in the caption of figure 01 
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Figure 6 : Diffractive dijet ep cross sections in the photoproduction kinematic range differential 
in zip, xjp, y and x 7 . The inner error bars represent the statistical errors. The outer error bars 
indicate the statistical and systematic errors added in quadrature. The overall normalisation 
uncertainty of 6 % is not shown. NLO QCD predictions based on the H12006 Fit-B DPDF set 
and the GRV 7 -PDF set, corrected to the level of stable hadrons, are shown as a white line. 
They are scaled by a factor 0.83 to account for contributions from proton-dissociation which 
are present in the DPDF fit but not in the data. The inner, light shaded band indicates the size 
of the DPDF uncertainties and hadronisation corrections added in quadrature. The outer, dark 
shaded band indicates the total NLO uncertainty, also including scale variations by a factor of 
0.5 to 2. A variant of the NLO calculation using the AFG 7 -PDF set is shown as a dashed line. 
For each variable, the cross section is shown in the upper panel, whereas the ratio to the NLO 
prediction is shown in the lower panel. 
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Figure 7: Diffractive dijet ep cross sections in the photoproduction kinematic range differential 
in E^ ctl , M x , A// jetK and (// Jef " s ). The inner error bars represent the statistical errors. The outer 
error bars indicate the statistical and systematic errors added in quadrature. Further details are 
given in the caption of figure [6] 
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Figure 8 : Diffractive dijet cross sections in the 7 p- and in the DIS regime normalised to the 
NLO calculation as a function of the photon virtuality Q 2 . The inner error bars represent the 
statistical errors. The outer error bars indicate the statistical and systematic errors added in 
quadrature. The data points are displayed at the geometrical bin centre. The NLO QCD pre¬ 
dictions are based on the HI2006 Fit-B DPDF set and, in case of photoproduction, on the GRV 
7 -PDF set, corrected to the level of stable hadrons. They are scaled by a factor 0.83 to account 
for contributions from proton-dissociation which are present in the DPDF fit but not in the data. 
The inner, light shaded band indicates the size of the DPDF uncertainties and hadronisation cor¬ 
rections added in quadrature. The outer, dark shaded band indicates the total NLO uncertainty, 
also including scale variations by a factor of 0.5 to 2. Also shown is the ratio of the RAPGAP 
MC to the NLO prediction. 


32 





















w HI VFPS data 


NLO HI 2006 Fit-Bx 0.83 x (1+5 ) 

hadr 



0.5 1 1.5 

ratio to NLO 


Figure 9: Diffractive dijet DIS and photoproduction cross sections normalised to the NLO cal¬ 
culation. Also shown is the double ratio of photoproduction to DIS cross sections, normalised 
to the corresponding ratio of NLO predictions. The inner error bars represent the statistical er¬ 
rors. The outer error bars indicate the statistical and systematic errors added in quadrature. The 
NLO QCD predictions are based on the HI2006 Fit-B DPDF set and GRV 7 -PDF, corrected 
to the level of stable hadrons. They are scaled by a factor 0.83 to account for contributions 
from proton-dissociation which are present in the DPDF fit but not in the data. The inner, light 
shaded band indicates the size of the DPDF uncertainties and hadronisation corrections added 
in quadrature. The outer, dark shaded band indicates the total NLO uncertainty, also including 
scale variations by a factor of 0.5 to 2. Variants of the NLO calculation, normalised to the de¬ 
fault NLO prediction, are also shown: the effect of using the AFG 7 -PDF parametrisation is 
studied in photoproduction. An alternative functional form of the scale is studied both in DIS 
and in photoproduction. 
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Figure 10: Ratios of diffractive dijet photoproduction to DIS cross sections differential in | A?/| 
and y. The inner error bars represent the statistical errors. The outer error bars indicate the 
statistical and systematic errors added in quadrature. The NLO QCD predictions are based on 
the HI2006 Fit-B DPDF set and GRV 7 -PDF, corrected to the level of stable hadrons. They are 
scaled by a factor 0.83 to account for contributions from proton-dissociation which are present 
in the DPDF fit but not in the data. The inner, light shaded band indicates the size of the 
DPDF uncertainties and hadronisation corrections added in quadrature. The outer, dark shaded 
band indicates the total NLO uncertainty, also including scale variations by a factor of 0.5 to 
2. Variants of the NLO calculation, normalised to the default calculation, are also shown. An 
alternative functional form of the scale is studied differential in |A?/|. The effect of using the 
AFG 7 -PDF parametrisation is studied differential in y. 
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Figure 11: Ratios of diffractive dijet photoproduction to DIS cross sections differential in zp 
and Ej! et ] . The inner error bars represent the statistical errors. The outer error bars indicate the 
statistical and systematic errors added in quadrature. The NLO QCD predictions are based on 
the HI2006 Fit-B DPDF set and GRV y-PDF, corrected to the level of stable hadrons. They are 
scaled by a factor 0.83 to account for contributions from proton-dissociation which are present 
in the DPDF fit but not in the data. The inner, light shaded band indicates the size of the DPDF 
uncertainties and hadronisation corrections added in quadrature. The outer, dark shaded band 
indicates the total NLO uncertainty, also including scale variations by a factor of 0.5 to 2. A 
variant of the NLO calculation using the AFG y-PDF is shown as a dashed line. 
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